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Abstract. Under the assumption that the variations of parameters of nature and the current 
acceleration of the universe are related and governed by the evolution of a single scalar field, 
we show how information can be obtained on the nature of dark energy from observational 
detection of (or constraints on) cosmological variations of the fine structure constant and 
the proton-to-electron mass ratio. We also comment on the current observational status, 
and on the prospects for improvements with future spectrographs such as ESPRESSO and 
CODEX. 
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1. Introduction 



We propose to probe the dynamics of the 
equation of state of matter and energy in the 
universe by using observations of cosmolog- 
ically varying fundamental parameters. This 
requires a number of assumptions: first, that 
there is clear evidence for dark energy that 
could be attributed to a rolling scalar field or 
quintessence; second, that there is a cosmo- 
logical variation of the fine structure constant 
of magnitude Aa/a < 1CT 5 , as suggested by 
Keck/HI RES high resolution q uasar absorption 
spec tra jMu rphy et all (l2004i I2003I) . but see 
also ISrianand et al.l (I2007I) ): and thirdly, that 
any variation of a arises from the evolution 



of the quintessence field (IDvali & Zaldarriagal 
(l2002l) : lchiba & Kohril (120021) 1. 

Concerning the first assumption, there is 
currently little observational evidence for dark 
energy to be more than a bare cosmological 
constant. However, if dark energy indeed re- 
sults from an evolving scalar field, it could be 
expected to couple to other forms of matter and 
lead to variations of masses and couplings (im- 
plying the second and third assumptions) un- 
less some unknown symmetry principle explic- 
itly forbids these couplings. 

We take the coupling between the scalar 
field and electromagnetism to be X<*f = 
-jBfi^F^F^ where the gauge kinetic func- 
tion Bf((f>) is linear, Bf(<f>) = 1 -^k(<P -<po) (and 
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k 2 = 87rG). This can be seen as the first term of 
a Taylor expansion, and should be a good ap- 
proximation if the field is slowly varying at low 
redshift. Then, the evolution of alpha is given 
by 



Aa _ a — ao 
a ao 



= (k(4> - 4>o) . 



We can also consider the variation of fi = 
m p /m e . In grand unified theories we expect a 
correlation between the v ariation of a and fi 
dCalmet & Fritzschl d2002l) ) given by Afi//j = 
R Aa/a, where R is a model-dependent numer- 
ical factor arising from correlated variations 
of Aqcd, the Yukawa couplings, the vacuum 
expectation value of the Higgs field and a it- 
self. Under simple assumptions we obtain R ~ 
-20, which is in severe tension with observa- 
tions that indicate a nontrivial varia tion of a at 
high redshift but null var i ation o f \x dKing et al] 
(2008); iThompson et alj d2009l) ') with equal or 
better precision. This simple exercise illus- 
trates the potential of cosmological observa- 
tions of quasar absorption lines and variation of 
fundamental parameters in discriminating par- 
ticle physics models. 

2. Reconstruction procedure 

In order to test our third assumption we 
may verify whether a given model of 
dark energy proposed in the literature can 
fit the data. This has been don e by a 
number of authors dCopeland et alJ d2 004): 
lAnchordoaui & Goldberg! d2003l) : iDent et all 
d2009l) : iBento & Felipe! 120091) 1 and it was 
found that many models can satisfy all the con- 
straints, though only in islands of the parameter 
space. 

An alternative approach is to extract the 
quintessence scalar potential from the obser- 
vational evolution of the data o n the variation 
of a. In IParkinson et al.l (2004), it is assumed 
that both the form of the gauge kinetic func- 
tion Bf{(p) and of the equation of state param- 
eter w^(z) are known. The authors parametrize 
these functions and fit the parameters by com- 
bining quasar and Snla data. 

Her e we discuss a slightly different ap- 
proach dNunes & Lid sey (2004l): lAvelino et alJ 



(2006); lAvelinol dl008)). We parametrize 
Bf(<f>) with a linear function as above and 
parametrize the evolution of a and/or fi with 
a polynomial g{N) where N = -ln(l + z). 
Essentially all we need is a functional form of 
(p'(N), then we integrate the equation 



(1) cr' = -(k4>') 2 (o- + a~ 3 ) . 



(2) 



where <x = p^/poQ^o- The solution cr(N) then 
provides the evolution of the equation of state 
parameter through 



(3) 



There are thus three steps of the reconstruction 
procedure required before we can apply these 
equations. 

2.1. Step 1: Obtaining the data sets 

The first step consists of obtaining data on the 
evolution of a and/or fi, either from existing 
observations, or from simulated data, for the 
purpose of forecasting the accuracy of recon- 
structions with proposed future spectrographs. 
We will consider the second approach. The 
data are generated from the numerical evo- 
lution of the quintessence field for a specific 
scalar potential V ((/>). We consider the normal 
distribution with mean Aa/a = (k(4> - <f>o) 
where Aa/a(z = 3) = -0.5 x 10~ 5 and we 
chose R = -6. We have assumed that with the 
ESPRESSO spectrograph for VLT, 200 systems 
will be found to determine a and 50 objects 
to determine fi with sensitivity 5 — 5 x 10" 7 . 
With the CODEX spectrograph for the E-ELT, 
we consider 500 systems measuring a and 100 
constraining fi with sensitivity 6 = 10~ 8 . 

2.2. Step 2: Fitting the data 

In previous works we chose to fit the data with 
a polynomial g(N) = Aa/a = g\N+g2N 2 + ...+ 
g m N"\ then the velocity of the field is simply 
up = g'/£. 

2.3. Step 3: Estimating £ 

The only missing ingredient is the value of 
f. We must estimate its value from indepen- 
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6 =5e-07, 8 =5e-07, 




Fig. 1. Reconstruction results for the scalar po- 
tential V((p) = Vo(e 10 ^ + e - 1 **). The dashed 
line represents the fiducial model, the solid line 
the best fit reconstruction and the dark and 
light bands the corresponding lcr and 2cr er- 
rors. Upper left panel: using a measurements 
with ESPRESSO; upper right panel: using a and 
yU with ESPRESSO; lower panels illustrate the 
same reconstructions with CODEX. 



dent observations such as Snla or weak lens- 
ing. More specifically, from the relation 



-1 + 



3Q d 



(4) 



and substituting for Kip' in terms of g'(N) we 
can obtain £ at any redshift. For example, at 
redshift z = we have 



3 Q^l + w ) 



(5) 



For typical values of Q^, () =s 0.7, wq ~ 
[-0.99, -0.6] and gl ~ 10 -5 we obtain f ~ 
10~ 7 - 10~ 4 which is comparable to bounds re- 
sulting from tests of the w ea k equi v alence prin- 
ciple dOlive & Pospelovl d2002l) : iDent et all 
(2009)). In Fig. Q] we illustrate a reconstruc- 
tion example for ESPRESSO and CODEX using 
a alone, and a and yi data in combination. 



3. The Rosenband bound 

A strong constraint on the current variation of 
a was obtained recently using atomic clocks 
dRosenband et alj d2008l) ) 



a/a = (-1.6 + 2.3) x 10~ w yr~ 



(6) 



This result rules out many models of 
quintessence with a monotonic evolu- 
tion of the field with a linear coupling if 
Aa/a{z — 3) ~ 10~ 5 , including the example 
of the previous section. There are of course 
a number of ways of evading these bounds. 
For instance, one may consider an oscillat- 
ing evolution: a potential with a minimum, 
V = Vo(exp(10K</>) + exp(-O.5K0)) would now 
satisfy the Rosenband bound because there is 
an oscillation of the field when this reaches the 
minimum of the potential. 

Alternatively one can modify the gauge 
kinetic function dMarra & Rosatil (2005)). For 
example if we consider the following gauge ki- 
netic function 



b p = i - a<p - w , 



(7) 



then our procedure would still apply but now 
there are two parameters which must be deter- 
mined using independent data. The following 
relations may be used: 



q g 

, 3 8 



w 



2g'\3(l+w) 



+ w(Q^-l) 



q 2 \g) 3£Vvv+l) 



(8) 
(9) 
(10) 



These require knowledge of the slope w'(z) and 
second derivatives of the polynomial, g"(z) at a 
given redshift. The reconstruction is therefore 
less accurate and specially difficult if wq « - 1 . 

4. What do the current data tell us? 

In this section we are going to be even more 
open minded by trying to dismiss our theo- 
retical prejudices and simply seek to take the 
current data at face value and understand what 
they might be telling us. 
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Fig. 2. Parametrization (fTTT i. Upper panel: 
comparison with binned data; lower left: evo- 
lution of the equation of state parameter w(z); 
lower right: the scalar potential. 




Fig. 3. Parametrization d!2l l. Upper panel: 
comparison with binned data; lower left : evo- 
lution of the equation of state parameter w(z); 
lower right panel: the scalar potential. 



Le t us take the unbinned Murphy et al. 
(2003) data which suggest a smaller value of a 
in the past. Considering the Rosenband bound 
and in addition the Oklo and meteorite analysis 
which put bounds of Aa/a = (0.7 ± 1.8)x 10~ 8 
at redshifts z = 0.14 dGould et all d2006l)) an d 
Aa/a _ (15+ 1.5 ) x 10" 6 (lOlive etaiT l2004): 
iDent et alj (120081) ) at z = 0.45, we may be 
compelled to consider a sharp transition in the 
value of Aa/a at about redshift z — 1 (see also 
iMortonson et alj d2009l) ). With this in mind we 
propose to keep the linear dependence of the 
gauge kinetic function Bf(<f>) but introduce the 
following parametrization for the evolution of 
the scalar field 



.(11) 



This parametrization corresponds to a field that 
evolves from a local maximum of the scalar po- 
tential, falls in a steep well and rises again ap- 
proaching another local maximum (see Fig. [2j». 
The velocity of the field is therefore decreasing 
today and a large vacuum energy is attained. 



An alternative parametrization is a function 
that allows for a variation of the field at early 
times 



N_ 



.(12) 



This second parametrization does not require 
that the field is initially at a local maximum, 
but instead allows an equation of state param- 
eter that approaches w(z) ~ at large red- 
shifts (see Fig. [3}. These forms of the poten- 
tial look distinctly unnatural, however we em- 
phasise that our objective here is to disregard 
theoretical prejudices and to use the data al- 
most blindly in order to uncover viable forms 
of the scalar potential. On the other had, this 
simple exercise highlights the importance of 
an independent observational confirmation of 
these variations. 

Performing a likelihood analysis using the 
first parametrization we obtain the contour 
plots shown in Fig. |4] for the amplitude of the 
transition in Aa/a, A = c£ and the width of the 
transition, A. We observe that by including ex- 
tra constraints such as the Rosenband bound, 
Oklo and meteorites, the contours are tighter. 
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Fig. 4. Likelihood analysis for parameters A 
and A using only quasar data (upper panel) 
and quasar, Oklo and meteorites data and the 
Rosenband bound (lower panel). The several 
shaded regions represent lcr, 2cr and 3<x confi- 
dence regions. 

Combining the data and constraints on the 
variation of a with Snla data we are then 
able to constrain c and therefore Indeed, to 
large values of f correspond small values of c 
and therefore the luminosity distance is to all 
effects indistinguishable form a ACDM sce- 
nario. The only upper limit on f comes from 
tests of the equivalence principle. Small values 
of £, however, give an evolution deviating sub- 
stantially from ACDM and one should be able 
to put lower bounds on this quantity with cos- 
mological data at redshift z > 1 , as illustrated 
in Fig.|5] 

5. Conclusions 

We have shown that under simple assumptions 
we can determine the nature of dark energy, 
not by fitting the parameters of a scalar po- 
tential to cosmological data, but by perform- 




0.5 1 1.5 





-F 









0.5 1 1.5 2 

z 

Fig. 5. The effect of different values of £ on the 
luminosity distance for a fixed value of A — 
c(. The dashed line represents a quintessence 
model with large £ which is completely de- 
generate with a ACDM cosmology and the 
solid line represents a quintessence model with 
small We used here the parametrization in 
Eq. (TQ. 



ing the inverse procedure, that consists in us- 
ing quasar data to reconstruct the potential. We 
have seen that the evolution of the equation 
of state parameter can in principle be found, 
subject of course to the precision of future 
data. These observations have profound impli- 
cations as the simple knowledge of the sign of 
w'(z) can help us to favour or discard freezing 
models of quintessence (increasing w(z) with 
increasing redshift), thawing models (decreas- 
ing w(z)) and k-essence models (typically also 
decreasing w(z)). 

This type of reconstruction directly probes 
the scalar field dynamics and may be carried 
out, with current data, to redshifts beyond z = 
4, far higher than the limiting value z = 1.7 
of Snla searches. Future astrophysical tech- 
niques may extend this to ev e n higher red- 
shifts dLevshakov et all (120071) : iKozlov et all 
(120081) 1. Moreover, the observations can be 
done from the ground and are consequently 
cheaper than satellite-based observations. 

Here, we have presented the reconstruction 
procedure for a minimally coupled scalar field, 
but other models with non-canonical kinetic 
terms, couplings to matter or multiple fields 
might have further interesting phenomenologi- 
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cal properties and therefore lead to alternative 
approaches. 
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